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Abstract 



We demonstrate that Monte-Carlo simulation is a practical tool to study nonperturbative 
aspects of supersymmetric quantum mechanics. As an example we study DO-brane quantum 
mechanics in the context of superstring theory. Numerical data nicely reproduce predictions 
from gravity side, including the coupling constant dependence of the string a' correction. This 
strongly suggests the duality to hold beyond the supergravity approximation. Although detail 
of the stringy correction cannot be obtained by state-of-the-art techniques in gravity side, in the 
matrix quantum mechanics we can obtain concrete values. Therefore the Monte-Carlo simulation 
combined with the duality provides a powerful tool to study the superstring theory. 
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1 Introduction 



Supersymmetric Yang-Mills (SYM) theories play prominent roles in theoretical particle physics. 
Especially, maximally supersymmetric theories are of crucial importance for superstring/M theory 
[H El O SI IS] - Given that most interesting questions can be answered only through nonpertur- 
bative study, it is important to construct theoretical frameworks for that. However, it is not a 
straightforward task because of the notorious difficulties of lattice supersymmetry (SUSY). 

Why is lattice SUSY difficult? In lattice gauge theory, gauge symmetry is kept unbroken. 
Therefore gauge symmetry breaking terms are not generated by radiative corrections. On the 
other hand, supersymmetry cannot be preserved completely, because the SUSY algebra contains 
infinitesimal translation, which is broken on lattice by construction. Therefore, even if a given 
lattice theory converges to a supersymmetric theory at tree level, SUSY breaking operators can 
be generated radiatively. In order to control the divergence one needs some exact symmetries 
at discretized level. In 4d M = 1 pure SYM, by keeping the chiral symmetry one can obtain 
the correct supersymmetric continuum limit [6p . In several extended SYM in less than four 
dimensions, by keeping a part of supersymmetries intact, SUSY breaking operators are forbidden 
at least to all order in perturbation theory [3J. (A similar statement holds for the Wess-Zumino 
model, and numerically tested in [9].) 

In one dimension (i.e. supersymmetric matrix quantum mechanics), the situation is much eas- 
ier. Because the theory is UV finite, one does not have to rely on exact symmetries and hence a 
simple momentum cutoff prescription works [10]. In fact, as demonstrated in [TO], the momentum 
cutoff method is much more powerful than a usual lattice regularization - the convergence to the 
continuum and restoration of the supersymmetry become much faster, and the Fourier accelera- 
tion, which reduces the critical slowing down, can be implemented without any additional cost. 
As a result it is possible to perform detailed Monte-Carlo study. In this talk we show the result 
of the Monte-Carlo simulation of maximally supersymmetric matrix quantum mechanics, which 
describes a system of multiple DO-branes, and compare it with dual string theory. Numerical 
data nicely reproduce not only predictions from supergravity but also the string a' correction. 
This strongly suggests the duality to hold beyond the supergravity approximation. In fact, in 
the gravity side, the stringy correction cannot be evaluated completely because of the lack of the 
knowledge on the higher derivative correction to the supergravity action. In the matrix quan- 
tum mechanics, however, we can obtain concrete values numerically. Therefore the Monte-Carlo 
simulation provides a new powerful tool to study the superstring theory. 

Although we concentrate on a specific model relevant for the gauge/gravity duality, the same 
method applies to other theories as well. Our message is Monte-Carlo simulation is a powerful 
and practical tool to study (SUSY-) quantum mechanics. 

This talk is organized as follows. In § [2] we briefly review the gauge/gravity duality. In 
§ [3] we explain the simulation method. Then in § H] we provide the simulation results and their 
interpretations. The materials treated here appeared previously in In order to avoid the 

repetition, we put emphasis on an explanation of the duality, so that the physical meaning of the 
simulation results become clearer. 

2 For recent numerical studies, see [7]. 
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2 The gauge/gravity duality in the DO-brane system 



The gauge/gravity duality conjecture [U[5] claims type II string theory on black p-brane background 
and maximally super symmetric super Yang-Mills theory inp+l dimensions are equivalent. When 
p = 3 this is known as the AdS^/CFT^ correspondence. In this paper we concentrate on p = 
case. The "derivation" of the correspondence is the same as p = 3 - both SYM and (weakly- 
coupled) superstring theory descriptions are valid in a certain parameter region of a coincident 
Dp-brane system, and hence they should be equivalent. U(N) super Yang-Mills is obtained in the 
near horizon limit of the iV-coincident Dp-brane system, where bulk degrees of freedom decouple. 
In the same limit the black p-brane solution in type II supergravity reduces to the near extremal 
solution, and at large- N and strong 't Hooft coupling the supergravity approximation is valid. 
The near extremal solution is given by [12} [T3] [5] 
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where U is radial coordinate perpendicular to the black brane, yii is the coordinate parallel to the 
brane, ^ls-p represents the spherical coordinate of the transverse directions and a constant d p is 
given by d p = 2 7_2p 7r^ 9_3p ^ 2 r((7 — p)/2). Uq represents the horizon of the black brane, which is 
related to the Hawking temperature T as 

T_(7-p)Ut P,,2 _, (3) 



Note that the gauge coupling constant gyu has dimension of (mass) ( - 3_p ^ 2 . As a result, an effec- 
tive coupling constant describing the black hole thermodynamics is a dimensionless combination 
T~( 3 ~ p ^\, where A = gy M N is the 't Hooft coupling. For this reason, when p < 3 "strong cou- 
pling" is equivalent to "low temperature". The dilaton (f>, and hence the string coupling g s = e^, 
depends on the radial coordinate, except for p = 3. In order for the supergravity approximation 
to be valid, the metric should be weakly curved so that a' correction is negligible, and g s must be 
small so that the closed string loop correction is small. Roughly speaking, the latter corresponds 
to the planar limit (N — > oo with fixed 't Hooft coupling) and the former is achieved at strong 't 
Hooft coupling (low temperature). 

On the contrary to the superstring, SYM is defined nonperturbatively. Therefore, if the 
gauge/gravity duality is correct, it provides a nonperturbative formulation of superstring theory. 
As we will see, this viewpoint is important not only conceptually, but also practically - once one 
puts SYM on computer, stringy correction to the supergravity, which is difficult to evaluate from 
the string theory, can be calculated numerically. We demonstrate it in § [3J 
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As an example of the correspondence, let us consider the ADM mass of the black brane per 
unit volume and the energy density of the gauge theory. They are predicted to be the same value, 
by identifying the Hawking temperature of the black brane to the temperature in the gauge theory. 
For p = 0, the prediction from the supergravity is 

FX -1/3 

= 7.4(TA- 1 /3)i4/5 ) (4) 

where the 't Hooft coupling A = Qy^N has dimension of (mass) 3 , T is the temperature of the 
system and E is the energy density. This relation is expected to hold at A r ~ 10 / 21 <C TA -1 / 3 <C 1, 
with which the supergravity approximation is valid around the horizon. Note that the black 
0-brane has a positive specific heat, in sharp contrast to the Schwarzschild black hole. 

When one applies AdS/CFT duality to solve strongly coupled field theory, the Gubser- 
Klebanov-Polyakov-Witten (GKPW) relation p3], which relates the anomalous dimension of the 
operators in CFT to mass of the supergravity modes, is very useful. Similar correspondence is 
proposed for p / 3 case as well. The basic idea of the GKPW relation is as follows. Let us con- 
sider a system of iV-coincident D-branes, which gives black p-brane background, and one probe 
D-brane put far away from others. In the supergravity, the effect of the probe is described as a 
perturbation around black brane background, specified by the boundary condition {hi}. In gauge 
theory side, it corresponds to U(N + 1) super Yang-Mills theory whose symmetry is broken to 
U(N) x U (1). The supergravity modes {hi} couple to operators Oi in the gauge theory. For p = 0, 
precise dictionary between field theory operators and supergravity modes can be determined by 
looking at how "Matrix theory currents" couple to bulk gravity modes [IS]: one considers the 
DO-brane action in weakly curved background, which is the Born-Infeld action, expand it around 
the flat space to linear order in supergravity modes and see how they couple to Matrix theory 
operators. For example, bulk metric couples to the energy-momentum tensor. This observation 
leads to p = analogue of the GKPW proposal [16] 

e -S SG [h] = / exp 



( I dtYsHmm) . (5) 



where Ssc[h] is the supergravity action in black zero-brane background as a functional of boundary 
values of the fields {hi}. The l.h.s. is evaluated by solving the classical equation of motion by 
imposing an appropriate boundary condition specified by {hi}. In order for this approximation 
to be valid, distance between two operators must lie in an appropriate range, so that only weakly 
curved and small g s region of the bulk contribute to the correlation function. From this one 
can calculate correlation functions. At zero temperature, operators corresponding to supergravity 
modes are expected to follow a power law [16] 



0{t)0(t')) cx |t _ , (6) 

although the theory is not conformant The region where the supergravity approximation is valid 
can be found as follows. As the separation between two operators becomes large, the supergravity 
mode propagates deeper into the bulk, so that it picks up the effect of large g s at small U. On the 



The power law appears due to the generalized conformal symmetry, which is explained below. 
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other hand, if separation is too small, contribution from large U region becomes important because 
the radial coordinate is related to the energy scale in the gauge theory. From these considerations, 
it turns out the supergravity approximation is justified at A -1 / 3 <C \t — t'\ <C \~ 1 / 3 N 10 / 21 [16j . 
Note however that this is just a sufficient condition and the result of the approximation may 
extends beyond this region. As we will see numerically, this is the case indeed. 

Another useful quantity which can be studied by using the duality is the supersymmetric 
Wilson loop [IT] [IS]. At finite temperature, the loop winding on temporal circle (Polyakov loop), 

W = Pexp(j dt(iA(t)+niXi(t))\ , (7) 

where n = (m, • • • , rtg) is an arbitrary unit vector, is an order parameter of the confinement- 
deconfinement transition [18] , Scalar fields appears because strings pull D-branes. In gravity side, 
(W) is identified to the exponential of the area of the minimal string world-sheet ending on the 
loop [UJ, 

log(W) = 1.89 (t\- 1/3 ) 3/5 + ••• , (8) 



where • • • are possible logarithmic and constant corrections. 

In AdS/CFT correspondence, CFT operators are put on the 'boundary' of the near horizon 
region, U — >• oo. This procedure seems to be subtle, however, because one has to take the near 
horizon limit in order to establish the correspondence and 'boundary' is the place where this 
procedure may fail. It is widely believed that the conformal symmetry saves the situation; if the 
correspondence holds at small distance, it extends to longer distance as long as the conformal 
symmetry exists both in gauge theory and in the supergravity. However, the system is not con- 
formal when p ^ 3. Then what can protect the correspondence? Actually the geometry (pQ) is 
invariant under the generalized scale transformation 

t,y\\ -> c~H,c-\, (9) 

U -> cU, (10) 
9ym -> c 3 ~ p g YM , (11) 

where c > 0. There is a counterpart in gauge theory, 

A^Xi -> cA^cXi, (12) 
t,x — > c~ 1 t,c~ 1 x, (13) 
9ym ~> (?- p g YM . (14) 

This generalized conformal symmetry serves as an alternative to the conformal symmetry [19j 
l2"m H6][2l]. Note that, for p ^ 3, gauge coupling is also scaled reflecting the fact that it is 
dimensionful. (The effective coupling AT - ( 3-p ) is invariant.) Therefore, in string theory, string 
coupling constant g s is rescaled. Hence it is not a "symmetry" in the usual sense; it should be 
interpreted as a transformation in the whole moduli of string/M theory. 
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3 Setup for numerical simulation 



The action of the maximally super symmetric matrix quantum mechanics is obtained formally by 
dimensionally reducing lOd super Yang-Mills theory to Id: 



S = ^j— f dtTr\-{D t Xi 
9ym Jo I 2 



1 



[Xi , Xj 



1 



(15) 



where Dt = dt — i [^4(i), • ] represents the covariant derivative with the gauge field A(t) being 
an N x N Hermitian matrix. It can be viewed as a one-dimensional U(N) gauge theory with 
adjoint matters. The bosonic matrices Xi(t) {i = 1, • • • ,9) come from spatial components of the 
lOd gauge field, while the fermionic matrices ip a (t) (a = 1, • • • , 16) come from a Majorana-Weyl 
spinor in lOd. The 16 x 16 matrices 7^ in (|15p act on spinor indices and satisfies the Euclidean 
Clifford algebra {74, 7j} = 2<5jj. When we study the energy density and the Wilson loop, we are 
interested in the finite temperature. Therefore, we impose periodic and anti-periodic boundary 
conditions on the bosons and fermions, respectively. The extent /3 in the Euclidean time direction 
then corresponds to the inverse temperature f3 = 1/T. For evaluation of correlators, we impose 
periodic boundary condition for both bosons and fermions. The coupling constant gy m m (|15p 
can always be scaled out by an appropriate rescaling of the matrices and the time coordinate t. 
We take gyM = (A = 1) without loss of generality. 

We take the static diagonal gauge A(t) = 4diag(ci!i, ■ • ■ , a^v), where a a can be chosen to 
satisfy the constraint max a (a a ) — min a (a a ) < 2tt by using the large gauge transformation with a 
non-zero winding number. We have to add to the action a term 



5 FP = -^21n 

a<b 



sin • 



a a - a>b 



(16) 



which appears from the Faddeev-Popov procedure, and the integration measure for a a is taken to 
be uniform. 

At finite temperature, we make a Fourier expansion 



A A-l/2 

r =-(A-l/2) 



(17) 



n=-A 



The indices n and r take integer and half-integer values, respectively, corresponding to the imposed 
boundary conditions^ Introducing a shorthand notation 



fciH hfcp=n 

we can write the action (fT5l) as 5 = 5b + 5f, where 



(19) 
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X^. n Xt-- A Tr([X u X 3 f 



When we impose periodic boundary condition for fermions, they are expanded as 



1 ab/,\ \ ^ lab itont 



(18) 
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■=-(A-l/2) L 



ru) — 



) 



(20) 



The fermionic action Sf may be written in the form S"f 



\M Aa r;Bfis^aT'4 , ps^ wh ere we have 



expanded ^} ar = Xm=i ^ar^ i n terms of U(N) generators t A . Integrating out the fermionic 
variables, one obtains the Pfaffian Pf.M, which is complex in general. However, we observe that 
it is actually real positive with high accuracy in the temperature regime studied in § 14.21 At 
very low temperature relevant for § 14.31 an d with periodic boundary condition adopted in § 14.41 
the phase of the Pfaffian does fluctuate violently. Here we simply neglect the phase and use 
|PfjM| = det(D 1//4 ), where V = M.^M., instead of PfAi. Surprisingly, it leads to the results 
expected by the dualitjU. 

To simulate this system treating fermions fully dynamically, we adopt the RHMC algorithm 
|22j . The trick of the RHMC is to use the rational approximation x _1//4 ~ &o + Sfc=i J+b k > w hich 
has sufficiently small relative error within a certain range required by the system to be simulated. 
(The real positive parameters and can be obtained by a code [23] based on the Remez 
algorithm.) Then the Pfaffian is replaced by |PfA1| = J dFdF* exp (— Spf)> where 



using the auxiliary complex bosonic variables F, which are called the pseudo- fermions. 

We apply the usual HMC algorithm to the whole system as described in |24| , except that now 
we introduce the momentum variables conjugate to the pseudo-fermions F as well as the bosonic 
matrices Xi and the gauge variables a a . 

When we solve the auxiliary classical Hamiltonian dynamics, low momentum modes tend to 
evolve larger amount compared to high momentum modes. Therefore, by taking the step sizes 
of the evolution of lower momentum modes larger, configuration space can be swept out more 
efficiently. This method is called the Fourier acceleration [25]. In the lattice gauge theory, in 
order to apply the Fourier acceleration one has to transform the configuration to the momentum 
representation. On the other hand, in the momentum cutoff method, because we are working 
directly in the momentum representation, the Fourier acceleration can be implemented without 
any additional cost. Thanks to this advantage the simulation cost can be reduced drastically. 

The main part of the computation comes from solving a linear system (T> + bk)x = F (k = 
1, • • • , Q). We solve the system for the smallest bk using the conjugate gradient method, which 
reduces the problem to the iterative multiplications of Ai to a pseudo- fermion field, each of which 
requires 0(A 2 iV 3 ) arithmetic operations if implemented carefully. The solution for larger 6/%'s can 
be obtained as by-products using the idea of the multi-mass Krylov solver [26j . This avoids the 
factor of Q increase of the computational effort. 

5 The agreement suggests the phase quenching is valid, but is not a "proof. To justify it one should look at the 
correlation between phase and values of physical observables [33] . 
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4 Results 



4.1 (In) stability of the black hole 

Because the theory is super symmetric, it has a flat direction along which scalar fields commute 
each other. As a result, partition function is divergent even at finite volume. Then how can we 
study this system numerically? 

String theory interpretation of the flat direction is obvious. Because scalar eigenvalues cor- 
responds to the positions of DO-branes, it corresponds to a gas of DO-branes. Once they spread 
widely, there is no force between them because of the supersymmetry and hence each DO-brane 
propagates freely. Black 0-brane, which is considered in the gauge/gravity duality, is a bound state 
of many DO-branes. In the matrix model, it is a bound state of scalar eigenvalues. Because the 
supergravity approximation is valid at large- iV, and there the black 0-brane is stable, we expect 
there is a stable bound state of eigenvalues at large-iV. At finite- iV, because the string coupling 
constant is finite, closed string can propagate and sometimes it escapes from the black 0-brane. 
In other words, scalar eigenvalues can run away from the bound state. Hence the bound state 
should be at most metastable at finite- N, i.e. black hole is unstable due to the quantum stringy 
effect. 

This is exactly what we observe in numerical simulation. For fixed N, we observe a metastable 
state at high and low temperature, while at intermediate temperature eigenvalues spreads quickly. 
To obtain expectation values, we use only configurations in metastable states. As N increases, 
the metastable state becomes stabler and it exists in wider parametric region. 

4.2 Energy vs ADM mass 

In order to calculate the energy density efficiently, we use a trick introduced in [24]. In the 
momentum cutoff prescription, one obtains [27] 



As we will see shortly, our data is precise enough so that we can evaluate the deviation from the 
supergravity limit quantitatively. In the dual gravity point of view, the deviation at large- ./V and 
at finite 't Hooft coupling corresponds to the string a' correction. Higher derivative corrections to 
the type IIA supergravity start with a' 3 order [28]. By a simple dimensional counting, it takes the 
form (a' /R 2 ) 5 , where R is the curvature radius of the black brane background, which translates to 
(TA" 1 / 3 ) 9 / 5 . Therefore, at large- N, correction to the supergravity limit E/N 2 = 7.41(TA _1 / 3 ) 14 / 5 
should be c(TA- 1 / 3 ) 23 / 5 , where c is an unknown constant. 

In Fig. [1] and Fig. [2] we plot the energy density obtained from our simulation [29j . Fig. Q] 
is a log-log plot of 7 AIT 14 / 5 - E/N 2 versus T. At sufficiently low temperature, a power law 
behavior can be seen clearly (i.e. data points form a straight line), which indicates that higher 
order corrections are negligible. (A deviation from a power law at very low temperature is merely 
a finite cutoff effect. Indeed by increasing the momentum cutoff from A = 6 to A = 8 data points 
go closer to the straight line.) By fitting the result with an ansatz 





= 7.41T 14 / 5 - cT p 



(23) 
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we obtain p = 4.58(3) and c = 5.55(7). If we fix p = 23/5 instead, we obtain c = 5.58(1). 
It strongly suggests that the duality is correct including the a' correction. Also, by assuming 
the duality is correct at this level, Monte-Carlo simulation provides very powerful tool to study 
the stringy correction to the black hole thermodynamics. Note that the coefficient c cannot be 
determined from gravity side, because the detail of the higher derivative correction to the type IIA 
supergravity is not known. Therefore our numerical result can be regarded as a prediction from 
gauge theory side. Usually people use the gauge/gravity duality to calculate something difficult 
(strongly coupled gauge theory) by dealing with easier problem (supergravity). However here we 
can use something difficult to study something more difficult (string theory). 




1 14 

Figure 1: The deviation of the internal energy jpE from the leading term 7.41 T s is plotted 
against the temperature in the log- log scale for A = 1. The solid line represents a fit to a straight 
line with the slope 23/5 predicted from the a' corrections on the gravity side. 



4.3 Wilson loop 

Next let us consider the Wilson loop [30] . It turned out that agreement with dual gravity prediction 
emerges only at much lower temperature compared to the case of the energy density. In such a 
low temperature region, because we have to take the cutoff A to be large, we could study only 
rather small values of N, up to A = 8. 

In Fig. [3l data points are shown in high and low temperature regions. At intermediate tem- 
perature, simulation with modest values of A is unstable because of the scalar instability. Note 
that we have evaluated (log \W\) instead of log(W) in order to reduce numerical error. We have 
taken the absolute value because at finite- A arbitrary phase factor can emerge due to the center 
symmetry and hence (W) = 0. At large-A, tunneling to different phase is suppressed and (W) 
agrees with (|W|) up to a fixed phase factor. That we have taken a logarithm before taking the 
expectation value can be justified at large-A, where fluctuation is suppressed. 

As can be seen from Fig. [31 logarithm of the Wilson loop behaves as 1.89T -3 / 5 + const. In 
fact it is difficult to distinguish constant and logT from our data; from stringy correction and 
quantum fluctuation of the string world-sheet, both constant and logarithmic corrections should 
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Figure 2: The internal energy jr^E is plotted against T for A = 1. The solid line represents the 
leading asymptotic behavior at small T predicted by the gauge-gravity duality. The dashed line 
represents a fit to the behavior (|23|) including the sub-leading term with p = 23/5 (fixed) and 
C = 5.58. 

arise. We consider the "constant" is actually sum of these corrections. 
4.4 Correlators 

Next we show simulation result for two-point functions [31]. Here we consider a class of operators 
given by 

1 

AT 

which couples to NS-NS 2-form and R-R 1-form. The exponent v predicted from the supergravity 
is [H] 

21 , . 

v = -. (25) 

In the simulation, first we obtain the correlation function in the momentum space, 

( J+(p) Jf(- P )) ( P = 27m//3; n = 0, 1, 2, • • • ). (26) 

Because of the existence of the momentum cutoff A, the result can be trusted only at small p. In 
order to estimate the correlators at large p, we assume the behavior 

JHpW(-p)) = ^ (27) 

where a is a constant, and we determine coefficients o by fitting a few points close to A. 

In order to obtain the correlators in the coordinate space, we perform the Fourier transforma- 
tion, 

(J?®J?{0)) = (j+(0)J+(0)) +^2co S (pt) (j+(p)J+(-p)) . (28) 

p>0 
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Figure 3: The plot of (log \ W\) for A = 1 against T 3 / 5 . The cutoff A is chosen as follows: 
A = 12 for JV = 4; A = 0.6/T for N = 6, 8; A = 4 for N = U; A = 6 for N = 17. The dashed 
line represents the results of the high-temperature expansion up to the next-leading order with 
extrapolations to N = oo, which are obtained by applying the method in [32 j . The solid line and 
the dotted line represent fits for A^ = 6 and Af = 4 respectively, to straight lines with the slope 
1.89 predicted from the gravity side at the leading order. 

We terminated the sum w.r.t. p at p = ^ x 1000. 

As a concrete example, we consider correlators at A^ = 3, A = 16, (3 = 4. For n > 12, we 
obtain the two-point function by extrapolating the data with the ansatz (|27p . For the fitting, 
n = 10, 11 and 12 are used. We evaluated the error by using the Jack-knife method, by dividing 
samples to four bins. The result after the Fourier transformation is shown in Fig. HI Straight lines 
represents a power law behavior proportional to \/t 2v+l , where v = 21/5 is the prediction from 
the supergravity. Surprisingly, the expected power appears at such a small value of N. In fact 
the power law is reproduced with remarkable precision, already at N = 2 [33]. Furthermore, the 
power law extends beyond the parameter region discussed in § [2j It is very interesting because it 
suggests the gravity prediction can shed a light on deep IR region, which is relevant for the matrix 
theory conjecture for M-theory pp. 

5 Conclusion and discussions 

In this talk we explained recent Monte Carlo simulations for maximally supersymmetric matrix 
quantum mechanics, in the context of the superstring theorjjf). The data reproduces predictions by 
the gauge/gravity duality precisely, and furthermore it provides a nontrivial prediction for stringy 
correction to black hole thermodynamics. 

There are many future directions. First of all, the same simulation techniques can be applied 
to other quantum mechanical system as well. It would be nice if one can learn nonperturbative as- 
pects of interesting theories which cannot be calculated analytically. It is also interesting to study 

6 The same system has been studied in [33] and qualitatively consistent results have been obtained. 
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Figure 4: Log-log plot of the correlator (J+(i) J+(0)) (i = 1, 2, 3, 4) for iV = 3, A = 16, = 4. 



theories in higher than one dimension. For example three-dimensional maximally supersymmetric 
Yang- Mills theory is realized [35] around fuzzy sphere background of the plane wave matrix model 
[36], and hence can be studied with the momentum cutoff methocQ- Another background in the 
same matrix model is argued to correspond to 4d M = 4 SYM in the planar limit [38, llj. In 
two dimensions, lattice formulation works without fine tuning [37} [39l |4"0] to all order in pertur- 
bation theory, and numerical studies so far confirms it works also at nonperturbative level. (For 
4-SUSY system, absence of fine tuning is confirmed in [41], by using dynamical fermion, for two 
independent lattice formulations. The conservation of the supercurrent has also been confirmed 
[42 j.) From string theory point of view, 2d SYM is related [43(. 144] to the black hole/black string 
transition [33]. Recent numerical study in this context can be found in [41 1 146] . Finally, 4d 
N = 4 SYM at finite- A 7 level is realized [17] by combining a matrix model technique [35] and two- 
dimensional lattic^E Together with [38], it provides nonperturbative tools to study AdSs/CFT^ 
correspondence. Monte-Carlo study of these models is of crucial importanc^]. 
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9 For another approach to AdSs / C FTn using simplified model see section V of [49] and references therein. 
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